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ABSTRACT 



To study the global morphology of the galaxy distribution in our local neighbourhood we calculate the Minkowski functionals for a 
sequence of volume limited samples within a sphere of 8 Mpc centred on our galaxy. The well known strong clustering of the galaxies 
and the dominance of voids and coherent structures on larger scales is clearly visible in the Minkowski functionals. The morphology 
of the galaxy distribution changes with the limiting absolute magnitude. The samples, encompassing the more luminous galaxies, 
show emptier voids and more pronounced coherent structures. Indeed there is a prominent peak in the luminosity function of isolated 
galaxies for Mb ~ -14, which at least partly explains these morphological changes. We compare with halo samples from a ACDM 
simulation. Special care was taken to reproduce the observed local neighbourhood as well as the observed luminosity function in 
these mock samples. All in all the mock samples render the global morphology of the galaxy distribution quite well. However the 
detailed morphological analysis reveals that real galaxies cluster stronger, the observed voids are emptier and the structures are more 
pronounced compared to the mock samples from the ACDM simulation. 

Key words, galaxy formation - cosmology - large-scale structure - galaxies - simulations 



1. Introduction 

We study the global morphology of the galaxy distribution in 
the local volume using Minkowski functionals and compare 
the geometry and topology of the galaxy distribution to halo 
samples from a ACDM-simulation. Only in our local neigh- 
bourhood we are able to observe objects down to very low 
absolute luminosities and only there we w ill be able to con- 
' struct an almost complete galaxy sample. iKarachentsev et af] 
I (1200 4') presented the essential version of such a catalogue 
. of neighbouring galaxies. They outlined the main structural 
elements of the local volume: the local sheet (part of the 
, local supercluster), the groups containing about 2/3 of the 
■ total local volume popul ation, and the local (TuUy) void 
' (and some other voids). Ti khonov & KarachentsevI (1200 6') and 
Tikhonov & Klvpin (2009) analysed these voids and the galaxy 
distribution within. To study the global morphology of the 
galaxy distribution we use the same data - an updated cata- 
logue of neighbouring galaxies (Karachentsev, private commu- 
nication). 

Cold dark matter simulations together with a cosmologi- 
cal constant (ACDM-simulations) are considered to match the 
distribution of M*-like galaxies quite well. However on small 
scales there are problems with the abundance of small mass 
aggregations. The ACDM model predict s thousands of dwarf 
dark ma tter halos in the local group (iKlvpin et al.L fT999aL 
iMoore e t al., 1999; Madau et al., 2QQ§ ) whil e only ~ 50 are ob- 
served. Recently ITikhonov & Klvpinr(l2009h found in ACDM- 
simulations a severe overabundance (a factor of 10) of halos in 
the voids compared to the observed number of galaxies in nearby 
voids. The overabundance of halos on small mass scales as well 
as in the voids is certainly a problem for the ACDM-model. 
In the present study, we investigate the spatial distribution of 
halos from a ACDM simulation compared to the galaxies ob- 



served in the local volume. For t he comparison with halo sa m- 
ples, we use the procedure from Tikhonov & Klvpin (2009) to 
map the circular velocities of the dark matter halos to luminosi- 
ties. This allows the construction of mock galaxy samples with 
the same number of objects and the same luminosity function as 
observed in the galaxy sample from the local volume. Hence, we 
assume that the ACDM overabundance can either be solved by 
the detection of new (very) low surface brightness galaxies, or 
by mecha nisms suppressing th e galax y formation in small dark 
halos (see ITikhonov & KlvpinI (|2009[) for a discussion and ref- 
erences). Our focus is on the global morphology of the galaxy 
distribution compared to the halo distribution. 



1.1. The galaxy sample 

Over the past few years searches for galaxies with distances 
less than 10 Mpc have been undertaken using numerous observa- 
tional data including searches for low surface brightness galax- 
ies, blind HI surveys, and N IR and HI observations of galax - 
ies in the zone of avoidance (IKarachentsev et all |2004 l2007h . 
The sample contains about 550 galaxies (see Fig.[l]l, the dis- 
tances to the galaxies are measured independent of the redshift, 
mostly with tip of the red giant brunch stars, cepheids, the Tully- 
Fisher relation, and some other secondary distance indicators. 
The galaxies in the local volume sample map the "real" 3D dis- 
tribution. The absence of the "fingers of God" effect simplifies 
the morphological analysis of the local volume. The distances 
have bee n measured with an accu racy of 8-10% for most of the 
galaxies (iKarachentsev et all 12004) . According to Karachentsev 
(2007, private communication) this sample is nearly volume lim- 
ited within an 8 Mpc radius from our galaxy, and is reasonably 
complete for galaxies with an absolute magnitudes Mb < -12. 
Some density gradients visible in the local volume sample are 
not necessarily caused by an incompleteness. They can be ex- 
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Fig. 1. The distribution of the galaxies in the local volume on the 
sky. 

Table 1. The galaxy samples considered. 



sample 


n.o. galaxies 


magnitude cut 


Iv8ml2 


315 


Mb < -12 


Iv8ml4 


205 


Mb < -14 


Iv8ml5 


133 


Mb < -15 


IvSmie 


96 


Mb < -16 



plained more naturally with structures in the local volume. The 
majority of groups resides in a thin layer with a distance of 
< 0.3 Mpc from the supergalactic plane. The filaments of the 
local supercluster fall mainly in the supergalactic plane, causing 
an excess of galaxies in this region compared to the outer regions 
perpendicular to the supergalactic plane. The overall complete - 
ness of this sample is discussed in iTikhonov & KlvpinI (|2009|) . 
They concluded that a significant incompleteness of the sample 
is improbable. There still exists the possibility that the obser- 
vational sample is incomplete, missing very low surface bright- 
ness galaxies free from gas, beyond the sensitivity of modern 
ground-based telescopes. The limiting magnitudes and the num- 
ber of galaxies in the subsamples considered are summarised in 
Table □ 

1.2. Minkowski functionals 

The genus of an isodensity surface is a well known method to de- 
scribe the topology of cosmological density fields ( Gott III et al.L 
Il987h . Minkowski functionals provide a unifying framework for 
the topology and the geometrical quantities like volume, surface 
area, and integrated mean curvature. They have been developed 
for the morphological characterization of t he large scale distri- 
bution of galaxies and galaxy clusters by 'Mecke et al.' ('19941) 
and have succes s fully been used in cosmology ( Kerscher et al], 
Il997t iKerscheii \200Ql), porous and disordered media, dewet- 
ting ph enomena, fluid and magneto hydrodynamics (see iMeck3 
(I2OOOI) for a review). 

In order to quantify the spatial distribution of the galaxies 
we decorate the galaxies with balls of varying radii (see Fig.|2l). 
Consider the union set - lJ;Io ^/-(x,) of balls of radius r 
around the galaxies at positions x,, i - 1,2, .. .N, thereby cre- 
ating connections between neighbouring balls. The global mor- 
phology of the union set of these balls changes with radius r, 
which is employed as a diagnostic parameter. It seems sensible 
to request that global geometrical and topological valuations of 
e.g. A,, are additive, invariant under rotations and translations. 




Fig. 2. Galaxies from the local volume decorated with balls of 
varying radii. 



and are continuous, at leas t for convex bodies. With these pre- 
requisites Hadwige^ ( Il957h could show that in three dimensions 
the four Minkowski functionals M^(A,), jj. - 0, 1,2, 3, give a 
complete morphological characterization of the body A^. The 
Minkowski functional Mo(A,.) simply its volume. Mi (A, ) is an 
eight of its surface area, M2(A,.) is its mean curvature divided 
by In^, and M^iAr) its Euler characteristics multiplied by ^. 
Volume and surface area are well known quantities. The in- 
tegral mean curvature and the Euler characteristic are defined 
as surface integrals over the mean and the Gaussian curvature 
respectively. This definition is only applicable for bodies with 
smooth boundaries. In our case we have additional contributions 
from the intersectio n lines and intersection points of the balls. 
iMecke et al.l (1 1994b discuss the extension for a union set of con- 
vex bodies. 

One may express Minkowski functionals in terms of n-point 
coiTelation functions, similar to a perturbative expansion. But 
the strength of an analysis with Minkowski functions is the di- 
rect quantification of the morphology. The additivity of the lo- 
cal contributions to the Minkowski functionals ensures the ro- 
bustness of the Mink owski functionals even for small point sets 
(iMecke et al.L Il994l) . In an analysis with Minkowski funtion- 
als, using up to 385 galaxy clusters from the Reflex sample, 
lKerscheretaT](l2001h were able to clearly distinguish the dis- 
tribution of observed galaxy clusters from a Gaussian point dis- 
tribution. As can bee seen from the Table [T] we are confronted 
with samples with a similar (small) number of points. Closely 
related to our analy sis is the invest i gation of the void sizes, 
as performed bv ITikhonov & Klvpini j2009l) for the same cata- 
logue of nearby galaxies. Indeed the void probability distribu- 
tion function is equal to one minus the volume density, the first 
Minkowski functional. Hence our analysis complements these 
investigation by additionally using the other three Minkowski 
functionals, to quantify the geometry, shape and topology of the 
galaxy sample. 
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Fig. 3. The volume densities m^(Ar) of the 
Minkowski functionals determined from the 
galaxy samples Iv8ml2 (black, solid), Iv8ml4 
(red, dotted), Iv8ml5 (blue, dashed), Iv8ml6 
(green, dashed dotted) For reference the m^(A;-) 
of a Poisson process with the same number 
density as Iv8inl2 are shown (one-cr area, thin 
black lines). 



2. Morphology of the local Volume 

The galaxy and halo samples considered are only complete 
within a sphere of 8 Mpc radius. We use boundary corrections 
as discussed in AppendixlAlto calculate the volume densities of 
the Minkowski functionals in an unbiased way. In AppendixlB] 
we show how Minkowski functionals describe the morphologi- 
cal features and briefly discuss some stochastic models. 

In Fig. [3] we show the volume densities ot^(A,) of the 
Minkowski functionals for a series of volume limited samples 
from the local galaxy distribution. For reference the functionals 
of randomly distributed points (Poisson process) are shown. 
Compared to randomly distributed points, the volume density 
OTo(A,) increase is considerably delayed for all of the galaxy 
samples. The empty space in between the strongly clustering 
galaxies fills up later than for the purely random distribution. 
The balls on the clustering galaxies already overlap for small 
radii, causing a slow rise and a reduced maximum of the surface 
density m\(Ar) x 8 compared to Poisson distributed points. For 
large radii the surface density of the galaxies is above the values 
for randomly distributed points. The galaxies cluster on low di- 
mensional structures and consequently the balls have more room 
to grow. In section 12.11 we will comment on the excess surface 
density in the samples with absolute magnitude B < -15, as its 
is visible for large radii. 

Again, due to the clustering, the density of the integral mean 
curvature m2(Ar)27T^ > increases slower and the maximum is 
reduced in comparison to a Poisson process. In a Poisson process 
we get completely enclosed voids leading to the strong negative 
signal. However from the galaxy samples only a week negative 
signal is seen. The concave structures are less prominent, as ex- 



pected for clustering on planar (super galactic plane) or even 
filamentary structures. 

For the radius r - Q the density of the Euler characteristic 
miiAo)^ equals the number density of the galaxies in the sam- 
ple. Hence the difference at r = is only a reflection of the differ- 
ent sampling. With increasing radius more an more balls overlap 
and the Euler characteristic decreases. Then tunnels through the 
structure are forming, giving a negative contributions to in^iA,). 
In the galaxy distribution only a small positive values for m3(A,) 
can be observer, strengthening the previous observation that the 
voids are not completely enclosed. 

2. 1 . Morphology changing with the absolute Magnitude 

Akeady from Fig.[3]we get the impression that the morphology 
of the galaxy distribution changes significantly if we include 
galaxies with absolute magnitudes Mg < -15 in our analysis. 
However we have to be more careful. Minkowski functionals 
calculated from points decorated with balls do depend on the 
number density of the point distribution. One can derive explicit 
expression in terms of high order correlation functions quan- 
tifyi ng this non-t rivial dependence on the number density (see 
e.g. Mecke'200(f). To compare the galaxy distributions with dif- 
ferent limiting magnitudes we generate samples with the same 
number of points. We randomly subsample the galaxy samples 
Iv8ml2 and Iv8ml4 to the same number density as in the sam- 
ple Iv8ml5. This allows us to compare the volume densities 
m^{A, ) for the galaxy sample with Mg < -12, Mb < -14 and 
Mb < -15 in Fig|4] The Minkowski functionals of the samples 
with Mb < -12 and Mb < -14 agree within the error bars. Also 
the Volume density and the Euler characteristic mostly agree 
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0.015 




Fig. 4. Volume densities m^(Ar) of the 
Minkowski functionals determined from the 
galaxy samples Iv8ml5 (blue solid) compared 
to the subsampled galaxy samples Iv8inl4 (red 
dashed) and Iv8inl2 (black, one-cr error from 
randomly subsampling). 



between all the samples, whereas the sample with Mg < -15 
shows an increased surface density mi{Ar) for radii from 2.2 
to 3.2 Mpc. The voids become significantly more emptier going 
from limiting magnitude Mb < -14 to Mg < -15. For the den- 
sity of the integral mean curvature m2{Ar) we see that the neg- 
ative contributions occur only for larger radii in the Mb < -15 
sample. These emptier voids are surrounded by concave struc- 
tures at larger radii. 



This behaviour can be explained by looking at the luminos- 
ity function of isolated galaxies in the samples. Fig.|5]shows the 
number of isolated galaxies in the local volume versus the ab- 
solute B magnitude. We use two selection criteria for isolated 
galaxies: first a geometrical one - an isolated galaxies may not 
have any neighbour within a distance of 0.3 Mpc, which is close 
to the virial radius of an M* gala xy. As a second c riterion we use 
the tidal index (Tl) of Karachen tsev et all (|2004|) . comparing a 
local dynamical timescale of the galaxy to the Hubble time. We 
consider galaxies with a tidal index less than zero as isolated. As 
can be seen from Fig. |5] the abundances of galaxies selected by 
either geometrical or physical criteria match. In both cases the 
same prominent feature shows up, a peak in the number of iso- 
lated galaxies with an absolute Magnitude Mb ~ -14. This peak 
explains that mainly isolated galaxies are lost in going from the 
sample Iv8ml4 with Mb < -14 to the more luminous sample 
Iv8ml5 with Mb < -15. It is reasonable to assume that iso- 
lated galaxies predominantly reside in voids. Then the voids are 
emptier in Iv8ml5, which is exactly what we observe with the 
Minkowski functionals. 




-20 -18 -16 -14 -12 -10 -8 




Fig. 5. The number of isolated galaxies shown for bins of ab- 
solute magnitude Mb- Either a geometrical neighbourhood cri- 
terion (0.3Mpc-isolated) or a dynamical criterion based on the 
tidal index (Tl) are used to determine isolated galaxies. 



2.2. Error estimates 

The position of a galaxy on the sky is known with a high ac- 
curacy. However, the radial distance is estimated with some dis- 
tance indicator To quantify the influence of the distance error on 
our estimates of the Minkowski functionals we randomise the ra- 
dial distance r using a Gaussian probability law with mean r and 
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standard deviation O.lr. We observe that the Minkowski func- 
tionals of the randomised samples follow closely the functionals 
of the galaxy sample, well within the error bars. We compare 
this measurement error to the statistical error of randomly dis- 
tributed points with the same number density. This Poisson error 
is approximately twice as large as the measurement error ob- 
tained from randomising the distances. Also the error from the 
subsampling (sometimes called Jackknife error, see Fig.|4]i has 
approximately the same amplitude as the Poisson error. 

3. Comparison with mock samples 

The local volume is certainly not a fair sample of the Universe. 
However, only within our local neighbourhood we are able to 
probe the galaxy distribution down to very low absolute luminos- 
ity. In the local volume we have prominent structural elements 
like the TuUy void and the supergalactic plane. The extraction of 
mock samples requires a careful selection of the position in the 
simulation box to find comparable features. Moreover the simu- 
lation has to be performed in a box large enough to allow those 
structures to form and, on the other side, with enough resolution 
on small scales to resolve the low mass halos. 

3.1. Simulation 

We use N-body simulation provi ded by A. Klypin do ne with the 
Adaptive Refinement Tree code (iKravtsov et aU[l997 ). The sim- 
ulation was performed in a box with side lengths 160/! ' Mpc 
for a spatially flat cosmological ACDM r nodel, which is com- 
patible with the 3rd yeai- WMAP data (ISpergel et al.L |2007|) 
with the parameters h = 0.73, 0,„ = 0.24, Qbar = 0.042, 
= 0.76, and the normalisation erg = 0.75. In this simula- 
tion the total number of particles is 1024^, the mass of a parti- 
cle is 3.18 • IO'^/i^'Mq, the spatial resolution is 1.2/i"'kpc, and 
the circular velocity of the smallest resolved halo is 27 km/s. 
We re-scaled all data (coordinates and masses of halos) to real 
units assuming //q = 72km/s/Mpc, which is close to the value 
reported from the WMAP analysis. Halo s are extracted usin g 
the bound density maximum halo finder (Klv pin et al.lll999bh . 
which detects bound halos with their subhalos and calculates 
several physical and geometrical parameters of these (sub)halos. 
We use the maximum circular velocity Vc to rank order the ha- 
los. The maximum circular velocity Vc can be determined more 
stable from the simulations and is easier to relate to observations 
as compared with the virial mass. We have a monotone relation 
between mass and circular velocity Vc - For reference, halos with 
Vc - 50 km/s have a virial mass of approx. 10'°Mq and halos 
with Vc = 20 km/s have a virial mass of approx. IO^Mq. 

3.2. Selection of a neighbourhood comparable to the local 
volume 

The reasonably big volume of the simulation box allows us to 
select mock s amples that mimic the loca l volume features more 
closely. As in lTikhonov & KlvpinI (l2009l) we use several criteria 
to select spheres with a radius of 8 Mpc from the simulation box. 

1. No halos with mass > 2 ■ 10^^ Mq reside inside the 8 Mpc 
sphere. 

2. The sphere must be centred on a halo with 1.5 • lO'^M© < 
mass < 3 ■ lO'^Mokm/s (the local group analog). 

3. The number density of halos with Vc > 100 km/s inside the 
8 Mpc sphere exceeds the mean value in the whole box by a 
factor ranging from 1.5 to 1.7. 



4. The number density of halos with Vc > 100 km/s found in- 
side a 4.5 Mpc sphere exceeds the mean value in the whole 
box by a factor greater than 4; 

5. There is no halo more massive than 5.0 • 10" with a dis- 
tance in the range from 1 to 3 Mpc). 

6. The central halos of other mock samples are more distant 
than 16 Mpc. There is no overlap between the samples. 

As discussed in lTikhonov & Klvpin (2009) samples selected by 
the above criteria may be considered as close geometrical and 
physical cousins of the real local volume. We could selected five 
such samples from the simulation box. Visual inspection shows 
that our mock samples look quite similar to the local volume 
galaxy sample. 

3.3. Mapping between luminosity and circular velocity of the 
halos 

When assigning luminosities to dark matter halos we follow 
the prescription of Conroy et al. (2006): First we build an or- 
dered list of all the galaxies in the observed local volume sample 
ranked by their luminosity. Then we rank order the halos in our 
local volume candidate samples by their circular velocity. Now 
the luminosity of the brightest galaxy is assigned to the halo with 
the largest circular velocity. Then the luminosity of the second 
brightest galaxy is assigned to the second biggest halo and so 
on. This procedure preserves the galaxy luminosity function and 
our mock samples have the same number density as the observed 
galaxy sample (see Table. [T]). By c onstruction we hav e a mono- 
tonic relation between Vc and MR. IConrov et al.l (12006) showed 
that this prescription reproduces also the clustering properties of 
the more massive galaxies in SDSS samples. Nonetheless this 
mapping is not perfect. Only halos with Vc > 30 km/s are left 
within the mock samples, whereas in the observed sample galax- 
ies wit h Mb ^ -12 may have a s ignificantly lower circular ve- 
locity dTikhonov & Klvpinil2009l) . However if we would assign 
higher luminosities also to halos with small circular velocities 
we would end up with too many halos. 

3.4. Morphological comparison 

In Fig.|6H8]we compare the morphology of the galaxy distribu- 
tion to the morphology of the corresponding halo samples from 
the ACDM simulation. We show the average and the one-o- 
range for the volume densities of the Minkowski functionals es- 
timated from the five mock samples. 

The surface area, the integral mean curvature and the Euler 
characteristic from the Mb < -14 galaxy sample are well re- 
produced by the corresponding mock samples as seen in Fig.|2] 
Only the volume density mo is slightly reduced on large scales, 
indicating emptier voids. 

This changes if we compare with the samples including less 
luminous galaxies with limiting absolute magnitude Mg < -12 
(see Fig.|6]l. The stronger clustering of the galaxies is visible 
in the steeper decrease of the Euler characteristic m^iAr) for 
small r and the reduced maximum of the integral mean curva- 
ture m2(Ar). The well known fact, that the voids are emptier in 
the real galaxy distribution, can be seen from the reduced mo(Ar) 
for large r. 

The more luminous halo samples with Mb < -16 also shows 
some morphological discrepancies compared to the correspond- 
ing galaxy sample (Fig. [8]). mo(A,) is reduced for large radii, and 
we conclude that the voids are emptier in the real galaxy dis- 
tribution. Also for larger radii, the halo samples do not show the 
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pronounced coherent structures visible in the galaxy distribution, 
as can be seen from the increased surface density mi{Ar) of the 
galaxy samples as compared to the halo samples. 

As we have seen in sect. 12. 11 the morphology of the observed 
galaxy distribution is depending on the limiting magnitude. To 
investigate this behaviour in the mock halo samples we proceed 
similarly to sect. 12. H and randomly subsample the halo samples 
corresponding to Mb < -12 and Mb < -14 to the same num- 
ber density as seen in the Mb < -15 sample. The morphology 
of the halo distribution only shows a week trend, still within the 
error bars, if we consistently select galaxies with higher circu- 
lar velocities. Hence, the physical mechanisms leading to these 
emptier voids for Mb < -15 in the observed galaxy distribu- 
tion (c.f. Fig.m are not captured by the ACDM model, and/or 
the mapping from the circular velocity to the luminosity is more 
complicated at the faint end of the luminosity function. 

4. Summary and Discussion 

To quantify the global morphology of the nearby galaxy dis- 
tribution we calculated Minkowski functionals from a series 
of volume-limited samples extracted from the Local Volume 
galaxy catalog ue (an updated versio n of the Catalogue of Nearby 
Galaxies from iKarachentsev et af] (|2004)). These samples are 
nearly complete for Mb < -12 within a distance of 8 Mpc from 
our galaxy. The strong clustering is dominating the morphology 
of the galaxy samples on small scales. On larger scales the empty 
voids and coherent structures show up in the Minkowski func- 
tionals. It seems quite obvious to attribute these morphological 
features to the Tully void, covering approximately a quarter of 
the sample, and the supergalactic plane. By extracting a series 
of volume limited samples with increasing luminosity we en- 
countered a clear change in the global morphology of the galaxy 
distribution. Going from a sample with Mb < -14 to the more 
luminous sample with Mb < -15 the voids become significantly 
emptier Indeed we saw a prominent peak in the luminosity func- 
tion of isolated galaxies with Mb ~ -14. We checked from the 
catalogue that all these galaxies are dwarf irregulars with signif- 
icant relative star formation rates and hydrogen mass. Excluding 
these galaxies from a sample makes voids significantly emptier 
It will be interesting to check if these galaxies form a unique 
population and to which extend selection effects or true physical 
effects cause the peak in the luminosity function. We checked 
that the errors in the distance indicators do not influence our re- 
sults. 

In the second part of the paper we compared the morphology 
of the galaxy distribution to mock samples from a ACDM sim- 
ulation. The position of our mock samples inside the simulation 
box was chosen to mimic the main features of our real galactic 
neighbourhood, and the luminosities have been assigned to the 
halos in such a way that the observational number density and 
luminosity function are preserved. The overall picture, with clus- 
tering on small scale, large voids and coherent structures on large 
scales can be seen also in the morphology of the mock samples. 
This agreement in principle allows us to look at the morphology 
in more detail. In spite of the careful selection of the mock sam- 
ples we still see morphological differences between mock and 
galaxy samples. The galaxies cluster stronger and the observed 
voids are emptier than in the mock samples. This observation is 
going beyond the well known overabundance of halos in ACDM 
models. By construction we have the same number of objects 
in the mock sample as in the galaxy sample, and we enforce 
the same luminosity function. Hence we see differences in the 
geometry and topology of the spatial distributions. Clearly, five 



moc k samples do not allow any elaborate statistical tests (see 
e.g. iBesag & Digglel[l977h . but we may calculate a rough es- 
timate of the fluctuations in the Minkowski functionals. These 
fluctuations are of the same order as observed for a Poisson pro- 
cess. Therefore we claim to see differences in the morphology, 
only if the Minkowski functionals are consistently outside the 
one-cr range. 

Certainly, our procedure of assigning luminosities to halos is 
oversimplified. Especially for low luminosities the mapping be- 
tween halo circular velocity and luminosity may not be straight- 
forward. With our mapping only halos with circular velocity 
Vc > 30km/s are left in the mock samples. This is in agree- 
ment with theoretical predic tions for the least massive halo th at 
can host a galaxy (see e.g. iHoeft et al.l (|2006|) : iLoebl (|2008|) ). 
However, in reality one observes galaxies with Mb in the range 
from -12 to -13 wit h significantly lower rotational velocities 
jTikhonov & Klvpinl |2009). If we would assign higher lumi- 
nosities to halos with small circular velocities Vc we either end 
up with too many mock galaxies, or we would have to drop some 
halos with high Vc to keep the number of mock galaxies equal to 
the observed number Such a procedure would break the mono- 
tonic relation between Vc and Mb- To generate larger emptier 
voids in the mock samples, one would have to additionally in- 
clude halos with small Vc residing in the vicinity of filaments 
and sheets, and drop halos mainly residing in voids. Also halos 
in voids with a small Vc should remain dim. A more intricate, en- 
vironment dependent modelling of the luminosities for the mock 
halos seems necessary. Especially if we also want to understand 
the dependence of the global morphology on the absolute mag- 
nitude cut in our galaxy samples. 

Another factor explaining the morphological discrepancies 
might be that all the five mock samples we could extract do 
not hold a void which is comparable in its size to the local 
(Tully) void. This allows two lines of arguments, either one 
blames the observations or the simulation. Looking at Fig.[T]we 
see that the zone of avoidance (our galaxy) cuts right through 
the Tully void. However, several targeted and blind searches of 
the zone of avoidance have been performed. It is unlikely that 
we actually miss enough galaxies to cut the void into halves 
jKarachentsev et aU 12004'). On the other hand the simulation 
was conducted in a box with a side length of I6O/1 ' Mpc. This 
size is considered as reasonably large, since we are only inter- 
ested in mock samples within a sphere of 8 Mpc radius. It re- 
mains an open question whether a ACDM simulation in a larger 
box could provide us with mock samples with big enough voids, 
especially when we consider the halo distribution at the low mass 
end of the mass function. 
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Appendix A: Calculating Minkowski functionals 



Kerscher et al 



used for the calculations of the Minkowski func- 
is an update d version of the code developed by 



Mecke et al 



The code 
tionals 

(Il997h . based on the methods outlined in 
1994 . The code is made available to the general 
public via http://www.math.lmu.de/ ~kerscher/software/ 

In the following we p resent the boundary corrections used by 
ISchmalzing et al.l (|1996[) . All galaxies in our samples are within 
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Fig. 6. Volume densities m^(Ar) of the 
Minkowski functionals determined from the 
galaxy samples with Mg < -12 (short dashed 
red lines) compared to averages over the 
corresponding halo samples (shaded Icr area). 
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Fig. 7. Volume densities m^(Ar) of the 
Minkowski functionals determined from the 
galaxy samples with Mb < -14 (short dashed 
red lines) compared to averages over the 
corresponding halo samples (shaded Icr area). 
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Fig. 8. Volume densities m^(Ar) of the 
Minkowski functionals determined from the 
galaxy samples with Mg < -16 (short dashed 
red lines) compared to averages over the 
corresponding halo samples (shaded Icr area). 




Fig. A.l. A two-dimensional sketch illustrating the geometry of 
the boundary corrections used. = 8 Mpc is the radius of the 
sphere Dr centred on our galaxy enclosing the local sample, r is 
the radius of the individual balls centred on each galaxy, and the 
black area is the set A,, n Z)^_r. 



a sphere of radius R - 8 Mpc around our galaxy, and similar for 
the halo samples. Let Dr be this spherical window containing 
galaxies. Ar - UZi ^i(0 is the union of balls B,(i) of radius r 
centred on the i-th galaxy, respectively. 

Considering only A,, we miss contributions from galaxies 
outside Dr, i.e. from galaxies not included in our sample. A well 



defined way is offered by looking at the intersection Ar n £)«_,■ 
(iMecke & Wagneiill99 1'). Using the shrunken window Dr^,. we 
make sure that all possible contributions are taken into account, 
as illustrated in Figure lA.ll In order to measure the boundary 
contributions we have to calculate the Minkowski functionals 
MfiiArHDR-r) of the intersection of the union of all balls with the 
window, and the Minkowski functionals M^(Ds_,.) of the win- 
dow Dr-i- itself. 

Now imagine the union set A,, of balls on all galaxies in a 
very large (but still finite) patch of the Universe, with Ar c A^. 
Integrating over all movements (translations and rotations) of the 
window within this large patch of the universe we obtain the 
spatial average 

(m^ (a, n = r M^[ArngDR^r)dg. (A.l) 

Here gDR^r is a translation and/or rotation of the shrunken win- 
dow, is the group of motions, and dg the suitably normalised 
Haar measure on this group. The integral over motions can 
be expressed as a sum using the principal kinematical formula 
(lHadwigeiilT957l) 




[Ar n gDR^r) = X U ^v(A,)M^_,(Z)«_,). (A.2) 

v={) ^ ' 



For a homogeneous and isotropic galaxy distribution the spa- 
tial average coincides with the expectation over different real- 
izations, and MftiAr n DR^r) is already an unbiased estimate of 
{m^j {Ar n DR-r)y This observation together with eq. ( |A.2t al- 
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Fig. B.l. The volume densities m^(A,.) of the Minkowski func- 
tional determined from a Poisson process (solid black lines) and 
points randomly distributed on random line segments (dashed 
red lines). Both models have the same number of points, and the 
same sample geometry as the galaxies in Iv8ml2. 



lows us to remove the boundary contribution of the window. For 
the volume density this is almost trivial: 



mo(A,.) 



Mo(A, n D«_,) 
Mo(Dr_,) 



(A.3) 



For all the Minkowski functionals estimates for the volume den- 
sities can b e obtained by applying the following recursive for- 
mulcfl jFava & Santalol [l978l) : 

M/A,nD«_,) M^_,(Dfi_,) 
niuiAr) = — > l'^\my(Ar)— . (A.4) 

This linear transformation of the observed Mf,{Ar n Z)R_r) facil- 
itates the direct comparison with the volume densities (a^) 
calculated from Point process models, namely the Poisson pro- 
cess. 



Appendix B: Morphology with IVIinkowski 
functionals 

In the following we give a qualitative discussion of the 
Minkowski functionals for a Poisson process decorated with 
balls of radius r and relate them to the observed morphologi- 
cal features. The Minkowski functionals of randomly distributed 
points (a Poisson process) are shown in Fig. IB. II 

With increasing radius r the volume is filled until reaching 
complete coverage where the volume density mo(A, ) reaches 
unity. The surface density mi (A, ) x 8 increases with the radius 
r, reaches a maximum, and finally approaches zero when all the 



we use the convention 2,',=, ^« = for j < i 



volume is filled up and no surface area is left. For small radii 
the balls grow outwards and the positive density of the integral 
mean curvature m2{Ar)2n^ > is a reflection of the mainly con- 
vex structure of the union set of balls. The volume starts to fill 
and the integral mean curvature reaches a maximum. For even 
larger radii we get negative contributions to the density of in- 
tegral mean curvature stemming from concave structures - the 
balls are growing inwards into the voids. Finally »J2(A,.) ap- 
proaches zero when no surface area is left. For the radius r = 
each ball (point) gives a contribution of one to the Euler char- 
acteristic and the density of the Euler characteristic m-i(Ao)^ 
equals the number density of the sample. With increasing ra- 
dius more an more balls overlap and the Euler characteristic de- 
creases. Then tunnels through the structure are forming, giving a 
negative contributions to m-i{Ar). As a clear signal of completely 
enclosed voids in the Poisson process, m3(Ar) shows a second 
positive maximum. 

Beyond the completely random model one may consider ge- 
ometrically inspired models, e.g. points randomly distributed on 
randomly placed line segments. Minkowski functionals for such 
a low dimensional Poisson processes embedded in three dimen- 
sions are discussed by Schmalzing & Diaferiol (|2000). As an ex- 
ample we show the Minkowski functionals for a model using in 
the mean 2.1 points, randomly distributed on randomly placed 
line segments with a length of 8 Mpc. As can be seen from 
Fig. IB. II some of the morphological features of the galaxy dis- 
tribution discussed in Section|2]can be reproduced. Specifically 
the stronger clustering and the suppression of completely en- 
closed voids is clearly seen in the Minkowski functionals. The 
choice of these parameters does not seem realistic. Choosing 
more points per line segment or using shorter segments leads to 
a largely different values for th e Minkowski functionals. Indeed 
ISchmalzing & Diaferiol (l2000 l) have shown that only with (ran- 
dom) mixtures of similar Poisson models for filament-, sheet- 
and field-galaxies one can reproduce the small scale morphol- 
ogy of the galaxy distribution. As can be seen from their Fig. Bl 
these mixture models are not appropriate on larger scales. This 
is expected, since randomly placed geometrical objects will not 
be able to model coherent large scale features. Other models 
are based on the hierarchy of correlation functions. The explicit 
incorporation of higher mome nts for the cons truction of point 
process models is discussed in lKersclierl (1200 ll) . At least for the 
lowest order models the Minkowski functionals can be given ex- 
plicitly (Kerscher et al., 2001). Models, motivate d from s tatisti- 
cal physics and spatia l statis tics are discussed by iMeckd (l2000h 
and ,Mecke & StovanI (l2005l) . 
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